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Phani K.V.V. Nukala 

Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6164, USA 



o 
o 

(N 

bj:): 

< 

(N 
o 

S 

> 



I 

o 



> 

O 

cn 

00 
O 

o 



X 



Stefano Zapperi 

CNR-INFM, SMC, Dipartimento di Fisica, Sapienza — Universitd di Roma, P.le A. Moro 2, 00185 Roma, Italy and 

ISI Foundation, Viale S. Severo 65, 10133 Torino, Italy 

Mikko J. Alava 

Laboratory of Physics, Helsinki University of Technology, FIN-02015 HUT, Finland 

Srdan Simunovic 

Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6164, USA 

We analyze the effect of disorder and notches on crack roughness in two dimensions. Our simu- 
lation results based on large system sizes and extensive statistical sampling indicate that the crack 
surface exhibits a universal local roughness of (loc ~ 0.71 and is independent of the initial notch size 
and disorder in breaking thresholds. The global roughness exponent scales as = 0.87 and is also 
independent of material disorder. Furthermore, we note that the statistical distribution of crack 
profile height fluctuations is also independent of material disorder and is described by a Gaussian 
distribution, albeit deviations are observed in the tails. 



I. INTRODUCTION 

The statistical properties of fracture in disordered me- 
dia are interesting for theoretical reasons and practical 
applications . An important theoretical issue is rep- 
resented by the scaling of crack surfaces. Experiments 
on several materials under different loading conditions 
have shown that the fracture surface is self-afhne Q and, 
in three dimensions, the out of plane roughness expo- 
nent displays a universal value of ^ ~ 0.8 irrespective 
of the material studied Q . The scaling regime is some- 
times quite impressive, spanning five decades in metallic 
alloys [4. In particular, experiments have been done in 



metals [5[, glass [6|, rocks [7[ and ceramics [SjJ covermg 
both ductile and brittle materials. Later on, a smaller 
exponent ( = 0.4 — 0.6 was observed at smaller length 
scales. It was conjectured that crack roughness displays 
a universal value of ^ ~ 0.8 only at larger scales and 
at higher crack speeds, whereas another roughness ex- 
ponent in the range of 0.4 — 0.6 is observed at smaller 
length scales under quasi-static or slow crack propaga- 
tion [4]. It was recently shown that the short-scale value 
is not present in silica glass, even when cracks move at 
extremely low velocities jOi]. In addition, in granite and 
sandstone, one only measures a value of 0.45 even at high 
velocities [10, ■ The current interpretation associates 
the value <^ ~ 0.8 with rupture processes occurring inside 
the fracture process zone (FPZ), where elastic interac- 
tions would be screened, and the value ( ~ 0.45 with 
large scale elastic fracture [1, [l3|- In two dimensions, the 
available experimental results, mainly obtained for pa- 
per samples, indicate a r oug hness exponent in the range 

c~ 0.6 -0.7 [il,[il,[il[I|. 

In this work, we investigate the influence of frac- 
ture process zone on crack roughness in two dimensions 



through two key variables: material disorder, expressed 
as a distribution in breaking thresholds, and pre-existing 
notches. Material disorder and the size of pre-existing 
notches play a significant role in determining the size 
of the FPZ ahead of the crack tip. When the disorder 
is weak, the size of the FPZ is small and the material 
fracture response is dictated by the stress concentrations 
around the notches. On the other hand, when the dis- 
order is strong, a relatively large fracture process zone 
is generated ahead of the crack tips. Similarly, the in- 
fluence of pre-existing notches on FPZ in the presence 
of disorder is non-trivial. This is especially the case 
when the initial notch size is small and disorder is suf- 
ficiently strong to allow for significant damage accumu- 
lation. As the damage starts evolving, multiple cracks 
develop, which in turn influence the stress concentration 
around the initial pre-existing notch. Even in the sim- 
plest case of non-interacting cracks, the stress flelds be- 
come additive and hence the proportionality with respect 
to inverse of square root of the initial notch size is lost. 
The presence of interacting cracks further complicates 
this scenario and the stress concentration around notches 
depends in a non-trivial fashion on the initial notch size. 
For large notches, the effect of disorder should be weaker 
since the fracture process is dominated by a single crack. 

The question we would like to address is how the rough- 
ness of the fracture surfaces depends on the material dis- 
order and the relative sizes of the pre-existing notches, 
given their influence on fracture process zone. Studies 
on random fuse model with uniform and power law dis- 
order have indicated that spatial correlations in the dam- 
age accumulated prior to the peak load (the maximum 
load before catastrophic failure) are negligible and that 
the damage is accumulated more or less uniformly up 
to the peak load [l3|- This suggests that the origin of 
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self-afBne roughness in the random fuse model should 
not depend on whether there is strong or weak disor- 
der since the spatial correlations are built in the system 
only at the final stage of macroscopic failure. Earlier 
studies that investigated the effect of disorder on crack 
roughness are controversial: based on two-dimensional 
disordered beam lattice simulations, Ref. [3 suggested 
a universal roughness exponent of ~ 0.86, whereas us- 
ing two-dimensional disordered fuse lattice simulations, 
Ref. [Ill argued against the universality of roughness 
exponent. Clearly, the situation warrants further inves- 
tigation especially in light of the role played by the FPZ 
in the current interpretation for different values of the 
roughness exponents. 



II. MODEL 

In this paper, we study the effect of disorder and 
notches on the crack roughness, by numerical simula- 
tions of the two-dimensional random fuse model (RFM) , 
where a lattice of fuses with random threshold are sub- 
ject to an increasing external voltage [1^ [2lj. The re- 
sults show that the roughness exponent does not depend 
on the breaking thresholds disorder strength and on the 
presence of a notch. We consider a triangular lattice of 
linear size L with a central notch of length oq. All of the 
lattice bonds have the same conductance, but the bond 
breaking thresholds, t, are randomly distributed based 
on a thresholds probability distribution, p{t). The burn- 
ing of a fuse occurs irreversibly, whenever the electrical 
current in the fuse exceeds the breaking threshold cur- 
rent value, t, of the fuse. Periodic boundary conditions 
are imposed in the horizontal directions {x direction) to 
simulate an infinite system and a constant voltage differ- 
ence, V , is applied between the top and the bottom of 
the lattice system bus bars. 

A power-law thresholds distribution p{t) is used by as- 
signing t = , where X G [0, 1] is a uniform random 
variable with density px{X) = 1 and D represents a 
quantitative measure of disorder. The larger D is, the 
stronger the disorder. This results in t values between 
and 1, with a cumulative distribution P{t) = t^^^ . The 
average breaking threshold is < t >= 1/{D + 1), and 
the probability that a fuse will have breaking thresh- 
old less than the average breaking threshold < t > is 
P(< t >) = + l)y/°. That is, the larger the D 

is, the smaller the average breaking threshold and the 
larger the probability that a randomly chosen bond will 
have breaking threshold smaller than the average break- 
ing threshold. 

Numerically, a unit voltage difference, V = 1, is set be- 
tween the bus bars (in the y direction) and the Kirchhoff 
equations are solved to determine the current flowing in 
each of the fuses. Subsequently, for each fuse j, the ra- 
tio between the current ij and the breaking threshold tj 
is evaluated, and the bond jc having the largest value, 
maxj^, is irreversibly removed (burnt). The current is 



redistributed instantaneously after a fuse is burnt imply- 
ing that the current relaxation in the lattice system is 
much faster than the breaking of a fuse. Each time a 
fuse is burnt, it is necessary to re-calculate the current 
redistribution in the lattice to determine the subsequent 
breaking of a bond. The process of breaking of a bond, 
one at a time, is repeated until the lattice system falls 
apart. 

Using the algorithm proposed in Ref. [l^l, we have 
performed numerical simulation of fracture up to sys- 
tem sizes L — 320. Our simulations cover an exten- 
sive parametric space of (L, D and ao) given by: L = 
{64,128,192,256,320}; D = {0.3,0.4,0.5,0.6,0.75,1.0}; 
and ao/L = {1/32,1/16,3/32,1/8,3/16,1/4,5/16,3/8}. 
A minimum of 200 realizations have been performed for 
each case, but for many cases 2000 realizations have been 
used to reduce the statistical error. 



III. CRACK ROUGHNESS 

Once the sample has failed, we identify the final crack, 
which typically displays dangling ends and overhangs 
(see Fig. [iJ. We remove them and obtain a single val- 
ued crack line hx, where the values of x G [0,L]- For 
self-afSne cracks, the local width, uu^l) = {J2x(^^ ~ 
(1/0 hx)"^)^^^, where the sums are restricted to re- 
gions of length / and the average is over different real- 
izations, scales as w{l) ^ K for I ^ L and saturates to 
a value W = w{L) ^ L'^ corresponding to the global 
width. The power spectrum S{k) = (hkh^k) / L, where 
hk = /la; expi(27ra;A:/L), decays as S(k) ~ 
When anomalous scaling is present [H, [2Jj UHl i the ex- 
ponent describing the system size dependence of the sur- 
face differs from the local exponent measured for a fixed 
system size L. In particular, the local width scales as 
w{£) ~ ^Cioc^^C-Cioc ^ gQ ^;jja,t the global roughness W scales 
as L'^ with ( > Qoc- Consequently, the power spectrum 
scales as S{k) - fc-(2Cioc+i)2,2(C-Cioc), 

In the following we investigate the influence of disorder 
D and initial notch size ao on crack roughness. Figure[2^ 
presents the scaling of local and global crack widths in 
systems with different disorder values and an initial rela- 
tive notch size of uq/L = 1/16. The slopes of the curves 
presented in Fig. [5^ suggest that the local roughness 
exponent Qoc = 0.71 and is independent of the disor- 
der. The global roughness exponent is estimated to be 
C = 0.87, and differs considerably from the local rough- 
ness exponent Qoc- The collapse of the data in Fig. ^jp 
clearly demonstrates that crack widths follow anomalous 
scaling law. The inset in Fig. [2j3 reports the data col- 
lapse of the power spectra based on anomalous scaling for 
different disorder values. This collapse of the data once 
again suggests that local roughness is independent of dis- 
order. A fit of the power law decay of the spectrum yields 
a local roughness exponent of (loc = 0.74. This result is 
in close agreement with the real space estimate and we 
can attribute the differences to the bias associated to the 
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FIG. 1: (Color online) A typical final crack in a system of size L = 320 with initial notch size ao = 80. Note that the crack 
shows dangling ends and overhangs, which are removed to obtain a single valued crack line. The initial central notch is not 
considered in the roughness calculations. 



methods employed |2a |. 

The influence of initial notch size on crack roughness 
is presented in Fig. [5^. The curves presented in Fig. [3^ 
represent the scaling of local and global widths for various 
notch sizes ap. Once again, the local roughness exponent 
is estimated to be Qoc = 0.71 and is independent of the 
initial notch size, whereas the global roughness exponent 
C = 0.87. Figure presents the data collapse of crack 
widths based on anomalous scaling law, which once again 
confirms that crack roughness follows anomalous scaling. 
The collapse of the power spectra in the inset of Fig. 
for different notch sizes confirms that the local roughness 
is independent of the initial notch size. A fit of the power 
law decay of the spectrum yields a local roughness expo- 
nent value of Qoc = 0.77. The close agreement of these 
results with the C,ioc = 0.72 obtained for the unnotched, 



strong disorder case [27[ indicates that the crack rough- 
ness is universal and is independent of disorder and initial 
notch size. The global crack width W however scales as 
W = {{{hb - hbf)Y''^ - with C = 0.87 ± 0.03, and is 
also independent of disorder and crack size. 

The scaling properties of the crack profiles h{x) can 
also be studied using the probability density distribution 
piAh{i)) of the height differences Ah{i) = [h{x + i) - 
h(x)] of the crack profile between any two points on the 
reference line (x-axis) separated by a distance £. Re- 
cently, there has been a debate over the scaling of this 
p{Ah{i)) distribution [H, HI, HO] , i.e., whether the scal- 
ing properties of p{Ah{i)) can be described by a single 
scaling exponent C or multiple scaling exponents are re- 
quired to describe the scaling of p{Ah{£)). The self-afhne 
property of the crack profiles implies that the probability 
density distribution p{Ah{£)) follows the relation 



p(A,.m)-.r-/(^) 



(1) 



whereas multiscaling of fracture surface implies that Eq. 
[l]is not valid. Although multiscaling oi p{Ah{£)) was ar- 
gued in Ref. [1^ , it has been shown in references [5^, f3(i| 
that p{Ah{£)) follows a self-affine monoscaling relation 
given by Eq. [T] and that multiscaling is an artefact that 



results at small scales due to the removal of crack pro- 
file overhangs. In the following, we investigate whether 
disorder has any infiuence on the scaling of p{Ah{£)). 

First, we present the probability distributions p(A/i(£)) 
of the height differences Ah{£) — [h{x+£) — h{x)] for vari- 
ous bin sizes € = 4, 8, 16, 32, 64 for a disorder oi D = 0.75, 
a system size of L = 320, and a relative crack size of 
ao/L = 1/16. Figure [3] shows the collapse of the central 
parts of the probability distributions p{Ah{£)) for larger 
£ values, but still within the local width scaling regime. 
The deviation for smaller £ values may be attributed to 
steps in the single- valued crack height profiles, which in- 
evitably arise due to the removal of overhangs on the 
crack surface. For larger £ values, the central parts of 
these distributions approach Gaussian, but clear devi- 
ations can be observed in the tails of the distribution 
from a Gaussian distribution. Similar deviations in the 
tails were observed for the uniform disorder case as well 

Second, we present the collapse of p{Ah{£)) distribu- 
tions in Figs. [3a)- (d) for various bin sizes £ and disor- 
ders. The collapse of these p{Ah{£)) distributions for var- 
ious disorders D at each bin size £ indicates that p{Ah{£)) 
distributions and the roughness exponent Qoc are unaf- 
fected by the material disorder. 



IV. DISCUSSION 

In summary, the evidence presented in this paper indi- 
cates that the crack surface roughness is unaffected by 
the material disorder and the presence of pre-existing 
notches. This can be inferred from the fact that ma- 
terial disorder, whether strong or weak, has a significant 
infiuence on the amount of damage accumulated prior to 
the peak load; however, the spatial correlations in the 
damage accumulated prior to the peak load are negligi- 
ble [13]. Indeed, Figs. [6l^a)-(d) show the snapshots of 
damage and crack propagation at peak load in typical 
fracture simulations of system size L = 320 having weak 
to strong disorder and small to large pre-existing notches. 
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FIG. 2: (Color online) (a) Scaling of local and global widths 
w{l) and W of the crack for different system sizes L = 
{64, 128, 256, 320}, disorder values D and a fixed ao/L = 1/16 
value (top). The local crack width exponent Cioc ~ 0.71 
is independent of disorder and differs considerably from the 
global crack width exponent — 0.87. (b) Collapse of 
the crack width data using the anomalous scaling law (bot- 
tom). Lc = {L — ao) is the effective length of the crack 
profile. Collapse of the data for a given disorder value im- 
plies that local and global roughness exponents are indepen- 
dent of disorder. The inset shows collapse of power spectrum 
S(k) using the anomalous scaling law with ^loc = 0.71 and 
— 0.87. The slope in the inset defines the local exponent 
via —{2(ioc + 1) = —2.48. (a)-(b) present a total of 20 data 
sets. 



FIG. 3: (Color online) (a) Scaling of local and global 
widths w(Z) and W of the crack for different sys- 
tem sizes L — {64,128,256,320}, notch sizes ao/L = 
{1/32,1/16,3/32,1/8,3/16,1/4,5/16,3/8}, and a constant 
disorder of D — 0.6 (top). Once again, the local crack 
width exponent (^loc = 0.71 is independent of notch size ao 
and differs considerably from the global crack width expo- 
nent ( — 0.87. (b) Collapse of the crack width data using 
the anomalous scaling law (bottom). Lc = {L — ao) is the 
effective length of the crack profile. The inset shows collapse 
of power spectrum S{k) using the anomalous scaling law with 
Ooc ~ 0.71 and = 0.87. The slope in the inset defines the 
local exponent via —{2Qac + 1) = —2.53. (a)-(b) present a 
total of 32 data sets. 



By following the damage growth process, one can easily 
see from these figures that there is very little crack ex- 
tension at peak load whether the material is strongly dis- 
ordered or weakly disordered. In addition, Fig. [S] shows 
that the FPZ can not be defined from a single damage 
snapshot, but it is necessary to average the damage over 
many realizations of the disorder. When this is done, we 
find that ^fpz depends strongly on disorder D (ranging 
from one to 12 lattice units for the values considered) 
but only weakly on ao and L [3l| . The independence of 



the roughness exponent on D suggests that self-affinity 
in the random fuse model is not related to the FPZ. 

This is further corroborated by our numerical simula- 
tions on a simplified random fuse model in which failure 
events form a connected crack thereby excluding damage 
nucleation in the bulk [2^. In this model, after break- 
ing the weakest fuse, successive failure events are only 
allowed on fuses that are connected to the crack. Oth- 
erwise, the rules of this simplified model strictly follow 
those of the usual RFM. Consequently, this model tracks 
only the connected crack along with its dangling ends in 




FIG. 4: (Color online) The logarithm of probability distri- 
butions p{Ah{£)) of the height differences Ah{£) = [h{x + 
I) — h{x)] of the crack profile h{x) for various bin sizes 
^ = 4, 8, 16, 32, 64. The results are for crack profiles obtained 
using a system size of L = 320, a disorder of D = 0.75 and 
a relative crack size of a^/L — 1/16. As a guide for the eye, 
we present a Gaussian fit for £ = 32. It can be seen that the 
central parts of the distributions are Gaussian, but the tails 
do not follow a Gaussian. 
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FIG. 5: (Color online) Collapse of the probability density 
distributions p{Ah{i)) of the crack profiles for various disorder 
values D = 0.3, 0.4, 0.5, 0.75. (a) £ = 32 (top left), (b) ^ = 64 
(top right), (c) £ = 96 (bottom left), and (d) £ = 128 (bottom 
right). The results are obtained for a system size of 1/ = 320. 



a disordered medium, and hence forms the most simpli- 
fied model to study the effect of disorder on crack rough- 
ness. As was shown in Ref. [29{ , this simplified model ex- 
hibits the same roughness and height-height correlation 
characteristics as that of conventional RFM. Even for 
this simplified model, the collapse of the power spectrum 
and local width of crack profiles for different disorders 
{D = {0.25, 0.4, 0.5, 0.6, 0.75, 1.0}) results in a roughness 
exponent value of Czoc = 0.7, which clearly demonstrates 





FIG. 6: (Color online) Snapshots of damage and crack evo- 
lution in disordered notched specimens of size L = 320. (a) 
D = 0.3 and ao = 20 (top left); (b) D = 0.6 and ao = 20 (top 
right); (c) D = 0.3 and ao = 80 (bottom left); (d) D = 0.6 
and ao = 80 (bottom right). 



once again that disorder is irrelevant for identifying crack 
roughness exponent. 

As also shown in Ref. [29|, when the branching of 
the cracks or damage within the fracture process zone is 
not allowed thereby limiting the crack extension to only 
the crack tips, a local roughness exponent of 0.5 is ob- 
tained. But, as soon as branching was present, the value 
of the roughness exponent was increased from Qoc = 0.5 
to Qoc — 0.7. This may also explain why there is no tran- 
sition in the value of the roughness exponent from strong 
to weak disorder in the presented simulations since dam- 
age is always present even for the lowest values of dis- 
orders considered. A similar conclusion was reached by 
Bouchbiner et al. studying a model for crack growth with 
damage nucleation [3^ . According to this work, as soon 
as a FPZ was introduced in the model the roughness ex- 
ponent increases from Qoc — 0.5 to a higher value. Thus 
it appears that the roughness exponent in two dimen- 
sions does not depend on the size of the FPZ, but only 
on the fact that a FPZ is present or not. In three dimen- 
sions the situation should be different since experiments 
suggest that the roughness exponent displays a crossover 
precisely at the FPZ size. It would be interesting to in- 
vestigate this issue using three dimensional simulations. 
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